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摘 要 : 针对 通信 带宽 受 限 的 网 络 环境 ， 引 入 一 种 基于 动态 事件 触发 机 制 的 数据 传输 策略 ， 研 究 了 一 类 非 线性 网 络 
化 系统 在 随机 网 络 攻击 下 的 有 限 频 域 故 障 检测 问题 。 首 先 ， 在 考虑 故障 灵敏 性 和 扰动 鲁 棒 性 的 前 提 下 ， 利 用 状态 增 
广 的 方法 将 原 系 统 的 故障 检测 问题 转换 成 昌 1H, 滤波 问题 ; 然后 , 在 考虑 该 区 有 界 非 线性 和 随机 网 络 攻击 的 情况 下 ， 
将 故障 的 有 限 频 域 特 性 考虑 到 已 性 能 指标 的 设计 中 ,并 结合 有 限 频 输入 特性 , 给 出 有 限 频 故障 输入 下 的 故障 检测 滤 
波 器 与 动态 事件 触发 机 制 的 联合 设计 算法 ; 最 后 ， 通 过 搅拌 俭 式 反映 器 系统 的 仿真 算 例 验证 了 该 方法 的 有 效 性 。 
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Dynamic event-triggered fault detection for networked systems in finite-frequency domain 
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Abstract: This paper investigated the fault detection problem in finite-frequency domain for nonlinear networked systems 
under stochastic cyber-attacks. To save limited network resource, this paper introduced a novel dynamic event-triggered 
scheme. Firstly, under the consideration of fault sensitivity and disturbance robustness, the paper converted the addressed fault 
detection problem into an auxiliary H_/H, filtering problem by augmenting the states of the original System and the fault 
detection filter. Taking sector bounded nonlinearity and stochastic cyber-attacks into consideration, the design of HH. 
performance index included the frequency characteristics of fault signals. Combined with finite-frequency input 
characteristics, the paper proposed the joint design algorithm for fault detection filter and dynamic event-triggered scheme 
under the finite-frequency fault input. Finally, a simulation example of stirred tank reactor system verified the effectiveness 
of the proposed method. 
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0 引言 间 非 线性 随机 系统 的 故障 检测 问题 ， 通 过 设计 自 适应 事件 触 

二 发 控制 率 ， 使 得 触发 函数 的 阔 值 随 着 测量 输出 的 变化 动态 调 
网 络 化 系统 是 指 被 控 对 象 和 传感器 、 控 制 器 、 执 行 器 以 ” 整 。 文 献 [17] 针对 离散 时 间 网 络 化 系统 ， 在 设计 故障 检测 观 
及 其 他 系统 组 件 之 间 通 过 共享 通信 网 络 连接 而 形成 的 复杂 控 测 器 的 基础 上 ， 构 建 了 基于 集 员 估计 的 残 差 评估 动态 阔 值 ， 
二 制 系统 。 随 着 通信 技术 和 计算 机 网 络 技术 的 高 速 发 展 ， 控 制 。 进一步 降低 网 络 资源 占用 且 满 足 更 为 灵活 的 系统 设计 要 求 。 
系统 也 向 网 络 化 与 智能 化 方向 转变 ， 网 络 化 系统 逐渐 被 广泛 。 此外， 在 网 络 信道 传输 数据 过 程 中 ， 系 统 可 能 会 遭受 网 络 攻 
地 应 用 于 实际 工业 领域 中 ， 如 移动 无 线 传感器 网 络 中 、 无 人  ” 击 。 与 周期 采样 机 制 相 比 ， 事 件 触发 机 制 下 按 需 发 送 的 数据 
机 中、 远程 医疗 中 等 。 然而 ,共享 通信 网 络 为 网 络 化 系统 带 来 。 包 数 量 减少 ， 但 对 系统 性 能 的 影响 更 大 ， 若 这 些 数据 包 被 网 
诸多 便利 的 同时 , 也 带 来 了 一 系列 挑战 , 如 网 络 时 延 由 、 数据 ，” 络 攻击 窜改 ， 则 会 破坏 系统 的 稳定 性 53。 因此， 在 针对 动态 
丢 包 辐 、 通 信 约 束 史 等 问题 ， 甚 至 引起 故障 。 近 年 来 ， 随 着 人 事件 触发 机 制 下 故障 检测 的 研究 中 设计 抵御 网 络 攻击 的 故障 


们 对 控制 系统 可 靠 性 与 安全 性 要 求 的 不 断 提 高 ， 针 对 非 理 想 ”检测 策略 至 关 重 要 。 

环境 下 网 络 化 系统 的 故障 检测 问题 已 成 为 控制 领域 的 热点 丰 需要 注意 的 是 ， 上 述 基 于 事件 触发 机 制 的 故障 检测 方法 

完 问题 [9 。 只 考虑 了 全 频 域 内 系统 的 整体 性 能 ， 然 而 实际 工程 中 的 故障 
在 男 一 个 研究 前 沿 ， 为 了 有 效 地 利用 系统 中 有 限 的 通信 ”信号 往往 发 生 在 有 限 频率 范围 内 ， 例 如 : 绥 变 故障 处 于 低频 


统 
资源 , 事件 触发 机 制 应 运 而 生疏。 其 主要 策略 是 通过 给 出 一 。” 域 ， 卡 死 型 故障 的 频率 为 零 趾 。 为 了 更 好 地 刻画 有 限 频 域 特 
个 预先 设 定 的 事件 触发 条 件 来 判断 是 否 将 测量 输出 传输 至 滤 ”性 ,文献 [19] 提出 的 广义 Kalman-Yakubovich-Popov (KYP) 
波 器 , 这 是 一 种 * 按 需 执 行 ” 的 非 等 周期 触发 方式 。 文 献 [11,12] ” 引 理 提供 了 一 种 处 理 故障 信号 频率 特性 的 有 效 方法 。 为 了 将 
在 故障 检测 方法 中 引入 了 静态 事件 触发 机 制 ， 当 最 新 触发 时 ”广义 KYP 引 理 的 适用 范围 推广 至 非 线性 系统 ， 文 献 [20] 而 
刻 的 测量 值 与 当前 测量 值 之 间 的 差 值 超出 给 定 闵 值 时 ， 更 新 究 了 基于 静态 事件 触发 机 制 的 有 限 频 域 故障 检测 问题 ， 并 将 
数据 传输 的 触发 时 刻 。 然 而 ， 上 述 静 态 事件 触发 机 制 的 低 触 。 非 线 性 误差 动态 系统 转换 成 一 类 线性 时 变 参 数 系统 来 处 理 。 
参数 灵活 度 限制 了 其 适用 范围 ， 因 此 有 学 者 在 此 基础 上 ， 然而 ， 目 前 基于 动态 事件 触发 机 制 的 非 线 性 网 络 化 系统 的 有 
提出 了 动态 事件 触发 机 制 [3"! 习 。 文献 [16] 考虑 了 一 类 连续 时 ” 限 频 域 故障 检测 问题 尚未 引起 关注 ， 更 不 用 说 考虑 网 络 攻 击 
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的 情况 ， 这 也 是 本 文 研究 的 主要 动机 和 出 发 点 。 

基于 以 上 讨论 ， 本 文 研究 了 随机 网 络 攻击 下 非 线 性 网 络 
化 系统 的 有 限 频 域 故 障 检测 问题 ， 主 要 贡献 体现 在 : a) 所 设 
计 的 故障 检测 滤波 器 同时 考虑 扰动 鲁 棒 性 能 和 五 故障 敏 
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30>0; 4 、?、m、m 和 = 均 为 给 定 的 非 负 实数 (bo 为 
当前 测量 输出 且 。%(D=Xb-y7 ， 表 示 上 一 触发 时 刻 的 测量 
偷 出 5K) 与 当前 测量 输出 yo 之 间 的 差 值 ， 

(kK)=y(k), kelk,kn), isN 。 只 有 当 y(o 满足 触发 条 件 (10) 时 ， 


感性 能 ; b) 为 了 进一步 节约 网 络 资 源 , 采用 一 种 具有 更 高 灵 
活 度 参 数 的 动态 事件 触发 机 制 ; c) 依据 有 限 频 域 不 等 式 的 时 
域 解释 PU， 直接 在 时 域 上 设计 有 限 频 域 故障 检测 滤波 器 。 


1 ”问题 描述 


1.1 系统 建 模 
考虑 离散 时 间 非 线性 网 络 化 系统 如 下 : 
dt 0 
y(K)=Cx(k) + Dw(K) 
其 中 ,x(K)eR" 为 系统 的 状态 向 量 , w(K)eR" 为 外 部 扰动 且 满 
足 [0,m) 范 数 有 界 ; f(K)eR” 为 待 检 测 的 有 限 频 域 故障 信号 
y(k)eR" 是 测量 输出 向 量 ， 系 统 参数 矩阵 和 4、B、F、C 和 
D 为 已 知 的 适 维 常数 矩阵 且 (C,4) 可 观测 。 非 线性 向 量 值 函 
数 ®(k,x()) 满足 如 下 扇形 有 界 条 件 : 
[Dk,x(k) — PK, KK) — Sx(K) -HET 
x [PK, x(K)) — DK, RK)) — Sx(k)— KK)]<O 
其 中 ，9 5, eR 为 已 知 常数 矩阵 且 @B(k,0)=0。 
故障 信号 所 处 的 频率 域 可 分 为 低频 段 . 中 频段 与 高 频段 ， 
统一 描述 为 如 下 集合 : 
Qu 人 {OeRwl0<O<0} (3) 
当 0<4-4<2z 时 ， 式 (3) 中 的 集合 ov 表示 中 频段 ，4 和 
42 分 别 为 中 频段 的 下 界 和 上 界 。 
当 4=-2=-4 时 , 式 (4) 中 的 集合 oo 表示 低频 段 ，29 为 低 
频段 的 上 界 。 


CO) 


Qa{0eRw||0|<0} (4) 
当 0=0,, =27-0, 时 ， 式 (5) 中 集合 Qo 表示 高 频段 ， 0, 
为 高 频段 的 下 界 。 


Qs{0e Rw||0l>0,} (5) 
定义 1 考虑 如 下 离散 时 间 系 统 : 
WK+D)= An(Kk) + Br(k) (6) 


若 在 输入 信号 v(k)eL,[0,%) 的 作用 下 ， 系 统 状态 


>w+D- (Pk + D) aT <(2sng 上] Su (#) (07) 


或 ex Dk Dem ni] +) ermn)] <0 (8) 
k=0 


或 Piattd -nltd -oT > (2sin] pL) (9) 


则 称 os4[0o) 分 别 为 对 应 于 [2, 9] 的 低频 输入 ， 
对 应 于 [4, 2] 的 中 频 输 入 或 对 应 于 (-o -0,]uv[9,, +o) 的 高 
频 输入 。 式 (7)~(9) 中 的 9 表示 故障 信号 fo 的 频率 ，% 和 
0 分别 表示 故障 频率 范围 的 边界 值 ， ,=(9i+0s)/2， 
0 =(0,, —06)/2 o 

注释 1 定义 1 借鉴 了 文献 21] 中 对 于 有 限 频 域 不 等 式 
的 时 域 解释 。 输 入 信号 的 有 限 频 域 特性 表现 为 对 系统 状态 的 
驱动 能 力 ， 如 : 频 域 处 于 [-0, 9] 的 低频 输入 (6 els[0,%) 满足 
不 等 式 (7)， 意 味 着 在 输入 "6) 的 作用 下 ， 系 统 状态 Wb) 的 相 
对 变化 速率 低 于 2sin(912) 。 
1.2 动态 事件 触发 机 制 的 建 模 

测量 输出 (6) 通过 带宽 有 限 的 通信 网 络 进行 传输 ， 为 了 
节省 网 络 资源 ， 采 用 如 下 动态 事件 触发 机 制 : 


kn = inf {k > kA9K) + oy (Ky(K) -eer (k) Qe, (Kk) < 0} 
Qk+D) = LK) + oyT (Ky(K) eeT(k)Qe, (Kk) 


Pc 


(10) 


7(k) 满足 : 


测量 输出 数据 才 会 被 传输 至 滤波 器 的 输入 端 。 图 1 给 出 了 动 
态 事件 触发 机 制 下 采样 和 释放 瞬间 的 时 间 序 列 ， 由 此 可 见 ， 
传感器 发 送 数 据 是 非 等 周期 的 。 


Oh lh 2h 3h 4h Sh 6h 7h 


传感器 节点 


执行 器 节点 


. 到 达 膀 间 
图 1 动态 事件 触发 机 制 下 采样 和 释放 瞬间 的 时 间 序 列 


Fig.1 The time sequences of the sampling and releasing instants 


〇 采样 瞬间 @ 释放 瞬间 


under the dynamic event-triggered scheme 
当 4 趋 于 0 时 ， 触 发 条 件 式 (10) 可 改写 为 
i inf {x >klowy (Ky(k) -eeT(k)Qe, (Kk) < 0} (11) 


注释 2 


可 知 式 (11) 退 化 为 静态 事件 触发 机 制 。 相 比 于 文献 [11,12] 
中 考虑 的 具有 固定 触发 阔 值 的 静态 事件 触发 机 制 ， 本 文 所 考 
虑 的 动态 事件 触发 机 制 更 具 一 般 性 与 灵活 性 。 
1.3 故障 检测 滤波 器 的 设计 

为 了 检测 系统 中 的 执行 器 故障 ,设计 如 下 形式 的 滤波 器 : 
£(k+D)= AC(K) + PEK) +L[IK) -CEK)] (12) 
r(k)=$(K)—CE(K) 
其 中 ， (kK) eR" 和 rsR" 为 滤波 器 的 状态 向 量 和 残 差 向 量 ， 
LeR”" 为 待 设计 的 滤波 器 增益 和 矩阵。 本文 考虑 网 络 中 存在 
虚假 数据 注入 攻击 ， 其 通过 帘 改 数据 包 内 容 而 破坏 数据 的 真 
实 性 和 完整 性 。 滤 波 器 端 在 随机 网 络 攻击 下 的 输入 信号 Xi) 
可 以 表示 为 


Fk)=[1 -Qa(K)ICg (x(K)) + ok) yk) (13) 
a(k) 为 服从 Bernoulli 分 布 的 白 噪 声 序列 。 


P{w( = 二 = 也 fx(o}= 区 
P{a(K)=0}=1-B{a(R)}=1-& (14) 
V{a(W)}=E{(a(t) -a))=a1-a) 


其 中 ， 2Z el0,1 是 已 知 常数 ，P{.} 表示 事件 发 生 概率 ，E{*)} 表 

示 数 学 期 望 ，V{*} 表示 方差 。 当 a(k)=0 时 ， 表 示 信 号 在 传输 

中 可 能 受到 网 络 攻 击 干扰 ， 其 概率 为 1-& 。 当 a(k)=1 时 ， 表 

示 信 号 在 传输 中 没有 遭受 网 络 攻击 ， 系 统 正常 工作 的 概率 为 

& 。Z& 越 小 网 络 攻击 出 现 的 概率 越 高 。 由 于 网 络 攻击 的 攻 

击 效力 是 非 线 性 的 ， 故 假设 g(x(o) 是 一 个 非 线性 函数 ， 且 满 

足下 列 二 次 约束 条 件 : 

8" (Kk,x(k))g(k, x(K)) < rx (kK)HTHx(K) (15) 

其 中 ，x >0 为 已 知 非 线性 函数 gCk,x(R)) 的 有 界 参数 ， 互 为 常 
数 矩 阵 。 

定义 E(K)=[xT (Kk) er(0oOT ~ e(k)=x(k) -E(k) AD(K)= PDK, x(K))— 

(Kk,X(1) ,并 结合 式 (1)(12) 和 (13), 通过 状态 增 广 可 以 得 到 如 

下 滤波 误差 系统 : 

E(kK+1)= (A -GK)A) ER +(B -GK)B, wk)+ 
Ff (P+(El -G(R)E,)e,(Kk) + ED(k, x(Kk)) + 
EAD(K) + (G+ G(K)G, )g(k, x(k)) 

r(k)=(Ci + G(R)C,) EK) +(D, +o)D,) w(K)+ 

(Es+G(K)Ec)e, (K+(Gs -dK)G, )g(k,x(k)) 
y(K)=CE(k) + Dw(k) 


Sl 


(16) 


四 


其 中 ， 上 表示 第 :个 触发 时 刻 且 态 =0 


9(k) 是 柔性 变量 


其 中 ，0 代表 适当 维 数 的 零 和 矩阵 ，Z 代 表 适 当 维 数 的 单位 矩 


阵 ， 其 余 参数 矩阵 为 


二 A 0 0 0 二 B 2 0 
4 = a ,A,= ,B= ,B, = 
(1-aLC A-LC ZC 0 B-aLD LD 


To ro [Dro rei, 0 
E=| |,E=| |,E= |,E= |,F= |,G=| 
We 


- [0 . i 
-Lc| G=[(&-DC C], C=[C 0], DP=2D , 户 =D ， 


E;=al ,Ec=1,G;=(-@C,G,=C, &(k)=a(k)-& 。 

本 节 的 目的 是 设计 动态 事件 触发 机 制 与 故障 检测 滤波 器 ， 
使 得 滤波 误差 系统 式 (16) 均 方 浙 近 稳定 且 满 足 H/F 性 能 指 
标 ， 即 满足 下 列 要 求 : 


a) 滤波 误差 系统 式 (16) 是 均 方 渐 近 稳定 的 。 
b) 当 joO=0 时 ， 系 统 满足 如 下 H, 性 能 指标 : 

SOF srw OF (17) 
并 使 得 7 尽 可 能 的 小 。 7 越 小 , 扰动 信号 对 残 差 越 鲁 棒 。 
c) 当 w(R)=0 时， 系统 满足 如 下 性 能 指标 

EOF PEO (18) 
并 使 得 8 尽 可 能 的 大 。Z 越 大 , 故障 信号 对 残 差 越 敏感 。 


1.4 残 差 评估 机 制 
结合 故障 检测 滤波 器 生成 的 残 差 ， 选 取 残 差 评 价 函 数 
7(o 与 阔 值 J 来 检测 故障 。 7(o 与 分别 为 


J(k)=E ero| (19) 
J = sup Jr 
Sb (20) 


其 中 : 7 表示 有 限 的 评估 时 间 长 度 。 
小 来 判断 故障 是 否 发 生 : 

J(k)> J 厨 故障 和 报警 
< J 一 无 故障 


通过 比较 7Go 与 mw 的 大 


QL) 


2 ”主要 结果 

本 节 以 线性 矩阵 不 等 式 的 形式 给 出 了 使 得 滤波 误差 系统 
式 (16) 满 足 ,性 能 指标 7 与 .性 能 指标 6 的 充分 条 件 , 并 通 
过 求解 凸 优 化 问题 得 到 滤波 器 的 参数 。 下 面 给 出 推导 过 程 中 


将 用 到 的 引 理 。 
引 理 12 对 于 一 个 向 量 值 函数 wxz(o) 和 已 知 常数 算 


阵 不 等 式 51,5, eR” ， 若 存在 : 
[®,x (KD) -Sx [BE xR) -Ssx(1)]<0 (22) 
则 有 如 下 不 等 式 成 立 : 
国 ,6 xX(k) 
[x?() Tk,x x a (23) 


其 中 ， $=(S7S,+S1S)/2，S$,=-(ST +S1)/2 。 
引 理 2P3 对 于 和 矩阵 4 ，Q=Q 7 和 P>0 ， 和 矩阵 不 等 式 
-CQ+47P4<0 成 立 当 且 仅 当 存 在 一 个 矩阵 G ， 满 足 : 
一 47ICG 

be po ml 

2.1 稳定 性 及 残 差 对 扰动 的 鲁 棒 性 条 件 
定理 1 给 定 满足 0<4<h、 加 +h%<1、0z<oz<o 和 >0 
的 实数 丸 、 如 、o,、o, 和 6 ， 如 果 存 在 正定 对 称 矩 了 泗 PB， 算 

阵 W 、X 和 滤波 器 参数 矩阵 工 ,使 得 如 下 矩阵 不 等 式 成 立 : 

[N,,]s; <0 (25) 
则 称 滤波 误差 系统 式 (16) 均 方 渐 近 稳定 ， 且 满足 H, 性 能 指 
标 ”。 其 中 ， 对 称 和 矩阵 Nj 中 的 非 零 项 如 下 : 
Ni=-Pu+rH™H-S , Ns=-P,, 
AU-EC)ICTK Ni =-A™Wy,—(— 5)CTE ,Nin =MGC™X ， 
Ns=CT No =-(4-DaCT , Nus=CCT , N,,= 


Q24) 


Ns =-5$, Ni = -A™W— 


Ni = QC™X 5 


» 


9 -Ps—8, 9 N27 = 人 9 
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Ns =-/A™Wi tuC™E, No= A™W+CIX ,Nya=C7, Nss=-77, 
Nso =—B™W -HABT 和 TAGCDT Nao=-BrW -BW+ EDIX  ， 
Ni =HACDIX ，Nio=CDIXE ，Nis=DT，Nu=aDr，Nas=CDT， 
Ni =-2e0 ，Nio=HUGX ，Nio=5E , Nin=udX ，Nip=6X ， 
Nu=z1I ， Nis=Gl ， Nss=(h+h -DI ，Ne=Ny=Nas=-T ， 
Noos=-Wy ， Nao=-W， ， No=-UiWy ， Nin=-W ， 
Nao =1(1 -2)CTX Naw =(1-2)CTT Nsn =-AUGCCTXK 

Naw=-GC™X ， Nus=(-aCT ， Nes=-6C™ ， Ns=Nin=Pn— 


-HT ， 
Naa=-(a+co) 7 ， 


证 明 


No = Ni =P, -Wa HWY, 

Ni = Nis1s =—I 。 

构造 如 下 Lyapunov-Krasovskii 泛 函 

Vi(k)=ET RPER) + AK) 
沿 系统 式 (16) 的 轨迹 求 偏 差 可 知 : 

E{AV,(R)}=E{V, (K+D}-E{V,(R)} = 

GK+DPe(K+D-e (RPe(K) + OK+D -HK) = 

[E00)(AT -GRAT)+w (RK) (Br -GK) BT)+ BK)+ 
eT(P) (Er+G(k)ET)+ g(x(k) (GT + GK)GT )]x 
P,[(A-&(K)A,) ER) +(B -GK)B,) w(K) + Bk) + 
(E+&(R)E,)e,(k) +(G +a(K)G,)g(x(K))]- 
GPER) + -DAR + oy (KyR) -ey (Ke, (kK) 


根据 扇形 有 界 条 件 式 (2) 和 引 理 1 可 得 : 
[Bk,x(K) = Sx [DK, xK) = Ssx(k)]= [DK, x(K) -SE ET 
[Bk,x 1) -SEé(K)] = (RAN) <0 


Nio0 = Naw = Ps-W—W™ 


(26) 


(27) 


(28) 


[A®B() — Sie(k)] [AGO — Sse(k)]= 


[A®O) — SEE0)] [AP(K) — SsEsé()]=7" (HAsn(k) <0 (29) 


其 中 ， w=[ET(K) wi(k) er) Bk) DTk Xk) ADTK) gT(x(k))], 


Gk)= 9092(k) ,$=(S7S, +ST9)12 §,=-(ST +S1)/2 。 


ErSE, 0 0 0 Es, 0 0 

* 000 0 00 

* *00 0 00 

A=| * **0 0 00 
米 米 米 六 米 0 0 

米 米 米 六 米 米 0 
[ESE, 0 0 0 Ers, 0] 

* 0000 0 0 
000 0 0 

As=| * **00 0 0 
米 米 米 米 0 0 0 

米 * 米 米 洲 I 0 

* * 米 米 * 米 0 

其 中 ， =[7 0],E=[0 1]。 


结合 式 (27)~(29)， 动 态 事件 触发 机 制式 (10) 和 引 理 2 可 推导 


Ee 
上 上 


E{A (= 
E{AV, (KR)} + g™(k,x(k)g (Kk, x(k)— gk, x(Kk))g (Kk, x(Kk)) < 
E{AV, (PD}+ rE (KHIHoE(R) 一 STOCK x(k)) 8 (kK, x(k))— 
7 (KA KR) (RAK) -Ee (KPe(K)+ 

(4+h -DIK + +0)y (Ky(K)— 28e; (Kk) Qe, (Kk) = 
TPE -S22b — diag{0,—7*1,0,0,0,0,0} -S33328 |y(K) 
当 w(k)=0 时 , 等 式 矿 (diag{0,-y?7,0,0,0,0ky(k)=0 成 立 。 根 
据 Schur 补 引 理 及 不 等 式 (30)， 可 知 1 一 S236 -S228 <0 


G30) 


成 立 的 充分 条 件 为 
| GD 
和 E33 0 
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=, 0 0 0 _E1$, -Er$, 0 | GC 0 E 000 6 C 0 E 000 6 
* _y21 0 0 0 0 0 aC, 0 dE 0 0 0 -66,| HacC, 0 dE, 0 0 0 -66G,|+ 
六 汪 _2e0 0 0 0 0 0 0 0 0 0 0 0 I 0 0 0 0 0 
二 中 ， =| * 兴 水 二 用 -DT 0 0 0 | ， 四 = 三 
其 中 ，S (H+h-D) 人 0 i i 
* 0 0 0 0 0 0 
六 * 米 洲 米 一 7 0 
* +* -220 0 0 0 0 
米 米 米 米 米 米 -I 
> - * #* (4+h -DI 0 0 0 |<0 (35) 
ee on I * 米 * * = 0 0 
全 4 C3 米 六 * * * 一 7 0 
Br -aB; DD, 米 米 米 米 米 米 _I 
Er -aEr 0 
BP ESE-BSE+eoHH ,So=| 0 0 |,sas=| 0 则 系统 式 (16) 满 足 有 限 频 域 #. 性 能 指标 目 olh >pI7Coh 。 
2 人 -1 0 0 
好 ~ 一 
Ei 冯 0 内 中 ， To=(a+o)CTC + rHIH, -HISH, ， Hs=|* -1 0 Sy 在 不 
LO7 G&G} L0 * * pl 
ee 同 频 域 率 下 的 取 值 如 表 1 所 示 。 
1 2 上 
Dr 表 1 集合 0 与 矩阵 在 不 同 频 域 的 取 值 
Er QEr Tab.1 0 And® fordifferent frequency ranges 
Su=|0 0 |, Ss»=diag{-Pr,-Pn)} , Ss=-(0+0) 7 ， 频 域 0 = 
0 0 [2 0 0 
低频 中 <O0 * 加 0 
3 0" | | | L -P -2co0s(0; )Q/ | 
[Gi -CCGT | Pp 0 er 
中 频 0<0<0, . 志 0 | 
Su=diag{-1,-1}, Ho=[H 0], &¢=Va(-a) 。 * * -Pp-2c0s(0,)0, 
故 有 3 -S23w36b -S332B -kdiag{0,-y°1,0,0,0，0,0}n(k) <0 a [号 -人 
a a 高 频 [0l=0, “总 0 
成 因此， 滤波 误差 系统 式 (16) 均 方 渐 近 稳定 。 [> * -p+2cos(0,)0,] 


为 了 证 明 滤波 误差 系统 式 (16) 在 零 初 始 条 件 下 满足 期 望 
的 五 ,性 能 指标 Y ， 考 虑 如 下 评估 指标 函数 : 


J_ = 卫 Er rt -yw Ow)| < 


k=0 


E Pr rpm yw) EVO + EfV, (oo)} = 
ey (32) 
E Er or -pwr wd + AV | 
k=0 
Nn KE Spb -SETb Su Sh Jy() 
根据 Schur 让 引 理 和 式 (31)， 可 得 S11 Sh — SB 一 
S42w3h <0 ， 则 评估 指标 函数 满足 六 <0 ， 即 
Fr) SY wr (RwR) G3) 
羽 此 ， 滤 波 误差 系统 式 (16) 满 足 期 望 的 忌 性 能 指标 。 
结合 Shur 补 引 理 和 引 理 2 可 知 ， 式 (31) 可 改写 为 如 下 算 
阵 不 等 式 ; 
[3 = 己 13 Cy 
>» 0 0 G4) 
* SB 0 
* 米 米 
L 44 
其 中 丰 -| WA W'B W'E 0 W'E, WE, ay 
六 ” ”LeEmrd -oaW'B, -eoW'E, 0 0 0 AW.'G, 
Ey =diag {P, —W—W',P,—W,—Wi'} 。 
| W W 
Nz a dl d2 二 11 12 x 
人 ' | | bn | Se 


可 得 式 (25) 成 立 。 定 理 证 毕 。 
2.2 残 差 对 有 限 频 域 故 障 的 敏感 性 条 件 

定理 2 给 定 满足 0<4<4h、4+4<1、0<a so 和 a>0 
的 实数 4、42 、a、m 和 = ,如 果 存 在 正定 对 称 和 矩阵 环 、2r ， 
使 得 如 下 不 等 式 成 立 : 


A F EE 0 EE E 6G A F EE 0 EE E GG, 
-64, 0 -bE 0 0 0 6G|3|-64, 0 -dE, 0 0 0 66, 
I 0 0 00 0 0 T 0 0 00 0 0 


证 明 首先 考虑 无 外 部 扰动 下 滤波 误差 系统 (16) 的 中 频 
段 和 情况。 分 别 对 式 (35) 左 乘 和 右 乘 矩阵 
(Kk) =[ET(R) Tk) eT(k) PT,x(k) APTR) gT(x(k) 37(D] 及 其 转 置 , 并 结 
合 迹 运算 关系 式 wQv=wr(Qvu")， 可 以 推导 出 : 

—é" (K)Pe(K)+E (K+DP E(kK+D-—r (Rr + PB fT KFRK) + 

tr[Q,; (enE(k+DETK) +e rm EK) ETk+1)— 

2cosO, EKET(K))] + rET KHIHER) — (36) 

8 (Kk,x(k)) gk, x(K) -7 (KA + A ) NK) + 

(1 + -Dok)+(o +o)y" (kK)y(K)—2ceT(k) Ye, (kK)<O 

考虑 到 非 线性 项 的 扇形 有 界 条 件 、 网 络 攻击 函数 的 二 次 约 
束 条 件 以 及 动态 事件 触发 机 制式 (2), 对 于 kelk,k4) ,可 推导 出 : 

ET(K+DP, E(kK+D) -eT(KP, SEK) + IK+D)— VK)— 
rr +P FR + [QO (er Ek DETK) + (37) 
ere(K)ET(k + —2c0s0, EK)ET KR))]<O 


选取 Lyapunov-Krasovskii 泛 了 水 为 VW(WD=eE (Pe(k) +9(k)， 
则 有 : 
Ar (D-rTOOr(O+D2FTCOOACOO+[O， (eeE(K+D 
-ET(K) te E(K)ET (k+l) -2c0s0, EK)ETK))]<O (38) 
在 零 初始 条 件 下 ， 对 式 (38) 两 边 同时 取 k 从 0 到 = 进行 
累加 可 得 : 


0 


Dr r+ BST + M) EO (39) 


eT] 


时 中 M:=V [-e re(k+DETK) -em EK) ET (K+) + 2cosO, E(k+ 


k=0 


DeT(k+D)] 。 


根据 欧 拉 公式 ，-M 可 以 在 计算 后 转换 为 定义 1 中 不 等 
式 (8) 的 左 侧 项 , 故 M 是 半 正 定 的 。 由 于 4 >0，, 当 式 (39) 成 并 
时 ， 式 (39) 左 侧 的 最 后 一 项 为 非 负 项 。 因 此 ， 可 以 得 到 
六 oz037wr， 即 定义 1 中 的 有 限 频 域 忆 性 能 指标 8 得 到 满足 。 

同 理 ， 令 9;=-0s=-0; 可 以 得 出 执行 器 故障 发 生 在 低频 
段 时 的 结果 ; 令 % =6, 和 04 =27x-9; 可 得 出 执行 器 故障 发 生 在 
低频 段 时 的 结果 ,其 证 明 过 程 与 上 述 中 频段 类 似 。 定 理 证 毕 。 
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注释 3 不 等 式 (35) 给 出 了 有 限 频 域 故障 敏感 性 的 设计 
条 件 , 是 广义 KYP 引 理 在 事件 触发 机 制 与 系统 存在 非 线性 项 博 
况 下 的 推广 。 
2.3 滤波 器 参数 计算 

结合 定理 1-2, 使 用 MATLAB 中 的 YALMIP 工具 箱 求解 
如 下 凸 优化 问题 ， 可 以 获得 最 优 滤波 器 参数 ， 以 及 所 1 


性 能 指标 yw 和 Zou : 


Pp, un -Bb 
Po > ci 
st. (25), (35) 
其 中 ， 滤 波 器 增益 矩阵 可 由 L=W-?X? 计算 求 得 。 
注释 4 为 了 降低 事件 触发 机 制式 (10) 的 数据 发 送 率 触 
发 阔 值 参数 ，s 应 取 较 小 值 ， 然 而 取 值 较 小 可 能 导致 定理 
1-2 中 的 约束 条 件 无 解 。 因此 , 为 了 在 保证 故障 检测 性 能 的 前 


提 下 获得 较 低 数据 发 送 率 ， 提 出 滤波 器 式 (12) 与 动态 事件 甬 


发 机 秆 1 式 (10) 的 联合 口 设计 算法 如 下 : 


算法 1 


b) 若 定理 1-2 中 的 线 
入 步骤 c); 否则 ， 重 


c) 计算 a=e+t+Ae， 


LL ， 重 复 步 又 c); 和 否则， 
d) 输出 动态 事件 触发 机 制式 (18) 的 参数 (2, 2) 、 


H_/H, 
e) 选 


性 能 指标 yn 


双 残 差 评价 函数 J(k) 与 


合 设计 
50 和 步 长 增 量 As 。 


性 和 矩阵 不 等 式 可 解 , 保存 6=&6、Q、 工 ， 进 
置 步骤 a) 中 参数 。 


由 


有 限 频 故障 输入 下 的 已 /六 故障 检测 滤波 器 与 
动态 事件 触发 机 制 的 于 
a) 初始 化 事件 触发 参数 4 、、a 、、 


若 定理 1-2 中 条 件 可 行 且 e>0， 更 新 5 、2、 


进入 步骤 d)。 


和 ss 。 
与 阀 值 /w ， 


并 通过 式 (21) 


决策 ， 判 断 故 障 是 否 发 生 。 
3 ”仿真 实验 与 结果 分 析 
本 节 以 搅拌 每 式 反映 器 系统 PC 为 例 来 说 明 所 提 方 法 的 


效 性 ， 人 


行 故障 检测 


[al a 

Ga —C)—kCae Rr 

dt 
dCs _E _E 

Ce thCae W haCae Wr (41) 

AB Be E 
aT _ -0- Ar eT) +kCsAHEC 忘 
dt C, DC 


其 中 , C a C4 为 进 料 浓度 ，Cs 为 生成 物 浓度 ， 
Cs 为 生成 物 浓度 ，7 为 反映 温度 ， 太 为 冷却 剂 温度 ,VV 为 反 


映 器 容积 AN 
能 ， 


器 系统 的 示意 图 , 选 


xD 表示 进 料 物 4 在 1 时 刻 的 浓度 ， 


F 为 容积 流量 ，2 为 液体 密度 ，E/R 为 反映 激活 
Ar 为 热 交 换 系 数 ，AH 为 反映 热 。 图 
取 x0D=[xa(@) ao wT 作为 状态 变量 ， 
2 表示 出 料 物 8 在 1 时 


刻 的 浓度 ， 


3(D) 表示 1 时 


刻 反 映 室 的 温度 。 图 


o 


反映 器 系统 的 示意 图 
取 步 长 At=0.1s x 


形 如 系统 式 (1) 的 状态 


拉 离 散 化 处 理 


寸 系统 式 (41) 进 行 欧 
空间 表达 式 : 


x(k+1)= Ax(k)+ Dk, x(Kk)) + Bw(k) + Ff (K) 


=Cx(k) 


其 中 ， 参 数 矩 阵 为 4 


-0.15 
F=| 007 | ,C=[L 0 


—0.01 


+ Dw(k) 


0.3872 0.0222 0.0183 0.0095 
=| 0.2444 0.3897 0.0007| , B=|0.6474|,， 


-0.0685 0.9711 0.4008 0.6779 


0],D=0.2 。 


仿真 过 程 中 , 假定 执行 器 故障 发 和 9 
=0. 
函数 @®B(x(K))=0.5(51 +5,)x(k) +0.5(S5, — 5 )x(k) 


取 参 数 4 =0.1， 


3, =01, 0.=2, &= 


-sin(x(k)) 


其 中 


， 可 得 


1 为 搅拌 等 式 反 映 


2 为 搅拌 等 式 


ChinaXiv 合 作 期 刊 
第 39 卷 第 9 期 


sin(x(k)):= diag{sin(xi(k)), sin(x,(k)), sin(xs(k))} , 9 ={0.1 0.15, 0.2} ， 
5S,={0.05, 0.1, 0.15} ， 五 =[1 0 0] ， 随机 网 络 攻击 的 函数 为 
g(x(KR))=—tanh[0.3x7(K) 0.2x7(Kk) 0.1x7(K)J",， 系统 的 状态 初 值 与 
滤波 器 初 值 分 别 为 x(0)=[0 -1 087 ，fzOD=[0 0 0F ， 外 部 扰 
动 信号 为 


w(k)=0.2e 0 sin(0.8k), 0<k<500 (42) 

通过 求解 凸 优化 问题 式 (40)， 可 以 得 到 最 优 性 能 指标 
/=00839 与 6=08194 ， 滤波 器 增益 和 矩阵 为 
工 =[-0.2529 0.0217 -04418 。 此 外 , 为 了 验证 本 文 所 提 方 法 的 
有 效 性 , 将 故障 检测 结果 与 不 考虑 故障 频率 特性 的 已 .滤波 器 
(参考 文献 [16]) 进 行 对 比 , 选取 与 上 述 有 限 频 域 方法 相同 的 7 


值 ， 即 y=00839 ， 求 得 滤波 器 增益 矩阵 为 Lm= 
[-0.4720 0.0217 -0.0795]? 。 
过 和 区 反应 基 吕 1 
pr 
三 | = 一 | 王 


搅拌 釜 式 反映 器 截面 图 
Fig.2 The cross-sectional diagram of stirred tank reactor 
本 节 主 要 考虑 两 种 类 型 的 执行 器 故障 ， 分 别 为 突变 故障 
和 绥 变 故障 。 首 先 ， 假 设 执行 器 在 时 刻 130<k<300 内 发 生 幅 
值 为 0.4 的 突变 故障 ， 表 达 式 为 
Lh 130<k<300 
0， ”其它 
图 3 给 出 了 残 差 信号 在 有 故障 和 无 故障 时 的 波动 情况 ， 
对 比 可 以 看 出 在 系统 发 生 故障 的 时 间 段 130< 上 <300 内 ， 残 差 
信号 有 较为 明显 的 波动 。 图 4 表示 在 突变 故障 发 生前 后 残 差 
评价 函数 /Co 随时 间 变 化 的 曲线 。 系 统 在 k=130 时 刻 发 生 故 
障 后 ，7(k) 逐渐 超过 闵 值 1%， 红 色 实 线 表 示 有 限 频 域 故障 检 
测 方法 下 的 残 差 评价 函数 ， 蓝 色 虚 线 表 示 无 故障 时 的 残 差 评 
价 函数 ， 根 据 残 差 评价 机 制式 (19)~(21) 进 行 计算 比较 可 得 
7434)=0.3649>J4=0.3643 , 即 4 个 时 间 步 长 内 检测 出 故障 的 发 
生 。 黑 色 点 划 线 表示 全 频 域 故障 检测 方法 下 的 残 差 评价 函数 ，6 


个 时 间 步 长 后 检测 出 故障 。 如 图 $ 所 示 , 在 500 个 评估 时 刻 内 ， 


图 2 


(43) 


一 由 


在 低频 域名 <z/110 , 选 
3, xk=1, 非 线性 


» 


采样 的 次 数 为 500 次 ， 而 传感器 发 送 数据 的 次 数 为 217 次 ， 意 
味 着 数据 传输 率 仅 有 43.4%， 节 约 了 56.6% 的 网 络 资源 。 
一 有 及 也 
一 一 -无 故障 J 


号 r(k) 


0.05 上 


残 差 信 


5 2 有 es 0 Te, 500 
时 间 (k) 
图 3 针对 突变 故障 的 残 差 信号 


Fig.3 Residual signal for abrupt faults 
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一 有 基隆 吕 下 的 有 限 大 城 广 法 
-- 有 歼 障 情况 下 的 全 频 域 方法 
了 情况 


四 


残 差 评估 函数 J(k) 


0 二 i 
时 间 (k) 


图 4 针对 突变 故障 的 残 差 评 


Fig.4 Residual evaluation function for abrupt faults 


价 函 数 


本 


50 


图 5 针对 突变 故 


100 150 200 250 300 350 400 450 500 


时 间 (k) 


障 的 触发 时 刻 与 触发 间隔 


Fig.5 Release instants and intervals for abrupt faults 


式 为 


0.005(k -130) 

f (Kk)=40.35+0.03sin(0.3x(k—200)), 200<k<300 
0， 

与 突变 故障 类 似 ， 
段 130<k<300 内 ， 


发 生 故 障 的 时 间 


的 波动 。 图 
时 间 变 化 的 


线 。 


区 
7 表示 在 组 变故 障 发 生前 
系统 在 上 =130 时 刻 发 生 故 障 后 ， 


色 实 线 表 示 


超过 阔 值 w ,经 


滤波 器 ) 下 的 残 差 评价 函数 ， 
居 残 差 评 价 机 财 


评价 函数 ， 根 ] 


J(135) =0.9547 > J,, = 0.9542 


其 次 ， 考 虑 执行 器 在 时 刻 上 =130 后 发 生 缓 变故 障 ， 表 达 


130 <k<200 
(44) 


已 


缓 变故 障 同样 处 于 低频 域 段 。 在 系统 
6 中 的 残 差 
站 后 残 差 评价 函 


言 号 有 明显 
数 /Co 随 
J(k) 逐渐 


有 限 频 域 故 障 检测 方法 (已 /已 


蓝 色 虚 线 表示 无 故障 时 的 残 差 
1(19)~(21) 进 行 计算 比较 可 得 
, 即 5 个 时 间 步 长 内 检测 出 故障 的 发 


生 。 黑 色 点 划 线 表示 全 频 域 故障 检测 方法 ( H, 滤波 器 ) 下 的 


残 差 评价 函数 ，12 个 时 间 步 长 后 检测 出 故障 。 
域 方法 19 采用 的 ,滤波 器 ,本 文中 
的 五 1 渡 ， 


可 知 ， 


更 快速 地 检测 


相 比 于 传统 的 全 频 
有 限 频 域 故障 检测 方法 选 
故障 的 频率 特性 ， 所 生成 的 残 差 对 故障 有 更 大 的 敏感 度 且 能 
出 故障 的 发 生 。 


结合 图 4 与 7 


器 能 够 充分 考虑 


残 差 信号 r(k) 


0 50 


100 150 200 250 300 350 400 450 500 


时 间 (k) 


图 6 ”针对 缓 变故 障 的 残 差 信号 


Fig.6 Residual signal for incipient faults 


器 发 送 数据 
节约 了 58.4% 的 网 络 资源 ,图 9、10 
E 事 件 触发 机 制 


周 


间 


残 差 评估 函数 J(k) 


此 外 ， 
8 所 示 ， 


期 性 


一 一 用 障 情 况 下 的 有 限 闫 域 方法 
4 一 一 有 故障 情况 下 的 全 频 域 方法 
-- - .无 故障 
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检测 


1 
300 


时 间 (k) 
图 7 针对 组 变故 障 的 残 差 评价 函数 


Fig.7 Residual evaluation function for incipient faults 


图 8 给 出 了 缓 变故 障 下 的 触发 时 刻 与 触发 间隔 。 如 

500 个 评估 时 刻 内 ， 采 样 的 次 数 为 500 次 ， 而 传 感 

的 次 数 为 208 次 ， 意 味 着 数据 传输 率 仅 有 41.6%， 
给 出 在 组 变故 障 的 影响 下 ， 


2 


个 触发 机 制 下 的 触发 时 刻 与 触发 


与 静态 事 


隔 。 


传输 议 


相应 4 
法 的 数 
文采 用 的 动态 事件 


中 


己 ， 


图 


Hi 


触发 时 刻 与 触发 间隔 


寸 比 医 


8~10 可 知 , 事件 触发 机 制 避 免 了 每 个 采样 周期 都 


8 


输出 数据 , 且 动 态 事件 


a 
© 


有 
已 


触发 机 制 能 进一步 降低 传输 率 。 
2 给 出 了 本 文 方法 与 不 同 传输 机 制 下 的 故障 检测 方 


传输 率 对 比 。 可 以 看 出 ， 在 突变 与 缓 变 故障 情况 下 本 


触发 机 制 均 能 更 显著 地 节约 网 络 资源 。 


| FLL 


0 50 100 150 200 250 300 350 400 450 500 
时 间 (k) 


事件 触发 机 制 的 触发 时 刻 与 触发 


本 文中 动态 间隔 


Fig.8 Release instants and intervals under the dynamic event- 


w 
a 


时 间 触 发 时 刻 与 触发 间隔 


图 9 文献 [9~10] 中 周期 性 时 


triggered scheme in this paper 
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a 


IN 


人 


ee ee 


Fig.9 Release instants and intervals under the periodic time-triggered 


故障 的 有 限 频 域 特 性 考虑 到 已 性 能 指标 的 设计 中 , 且 快 速 地 
检测 网 络 攻击 下 搅拌 釜 式 反映 器 系统 中 故障 的 发 生 。 
晰 地 呈现 了 本 文 仿 真实 验 数据 。 此 外 ， 引 入 的 动态 事件 触发 
传输 策略 可 以 进 
显著 地 市 约 网 络 资源 。 


更 


scheme in reference [9~10] 


综 上 所 述 ， 本 文 所 设计 的 五/ 故障 检测 滤波 器 可 以 将 


表 3 清 


步 地 减少 输出 测量 数据 向 滤波 器 的 发 送 量 ， 


录用 定稿 朱 洪 ， 等 : 


4 ”结束 语 


本 文 提 出 一 种 动态 事件 触发 机 制 下 非 线 性 网 络 系统 的 有 
限 频 域 故障 检测 方法 。 在 考虑 扇 区 有 界 非 线性 和 随机 网 络 攻 
击 的 情况 下 , 将 故障 的 有 限 频 域 特性 考虑 到 已 性 能 指标 的 设 
计 中 ， 利 用 Lyapunov 稳定 性 理论 和 线性 矩阵 不 等 式 方法 得 
到 故障 检测 滤波 器 存在 的 充分 条 件 ， 并 使 用 凸 优化 方法 获取 
最 优 滤波 器 参数 。 仿 真 结果 表明 ， 与 全 频 域 故障 检测 方法 相 
比 ， 有 限 频 域 方法 对 故障 更 加 敏感 ， 具 有 更 快 的 检测 速度 ， 
且 动 态 事件 触发 机 制 能 更 显著 地 节约 网 络 资源 。 此 外 ， 将 本 
文 所 提出 的 故障 检测 方法 推广 至 多 智能 体系 统 拘 与 信息 物 
理 系统 PC4， 是 下 一 步 需要 研究 的 问题 。 


人 


0 50 100 150 200 250 300 350 400 450 500 
时 间 (k) 


图 10 文献 [11~12] 中 静态 事件 触发 机 制 的 触发 时 刻 与 触发 间隔 


Fig. 10 Release instants and intervals under the static event-triggered 


Scheme in reference [11~12] 


表 2 不 同 传输 机 制 的 数据 传输 率 比 较 

Tab. 2 Comparison of different data transmission Schemes 
、 数据 传输 率 

ae 突变 故障 组 变故 障 
周期 性 时 间 触 发 机 制 100% 100% 
静态 事件 触发 机 制 " 六 58.1% 56.8% 
本 文 动态 事件 触发 机 制 43.4% 41.6% 

表 3 不 同方 法 下 故障 检测 性 能 的 比较 


Tab.3 Comparison ofthe fault detection performances by different methods 


故障 检测 时 间 


故障 检测 方法 ”五 , 性 能 指标 鼠 - 性 能 指标 


突变 故障 ”组 变故 障 
全 频 域 方 法 09 0.0839 一 6 个 时 间 步 长 12 个 时 间 步 长 
本 文 有 限 频 域 方法 “0.0839 0.8194 ”4 个 时 间 步 长 5 个 时 间 步 长 
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